UNCLASSIFIED 


Defense  Technical  Information  Center 
Compilation  Part  Notice 

ADPO 13709 

TITLE:  Shape-Measure  Method  for  Introducing  the  Nearly  Optimal 
Domain 

DISTRIBUTION:  Approved  for  public  release,  distribution  unlimited 


This  paper  is  part  of  the  following  report: 

TITLE:  Algorithms  For  Approximation  IV.  Proceedings  of  the  2001 
International  Symposium 

To  order  the  complete  compilation  report,  use:  ADA412833 

The  component  part  is  provided  here  to  allow  users  access  to  individually  authored  sections 
of  proceedings,  annals,  symposia,  etc.  However,  the  component  should  be  considered  within 
the  context  of  the  overall  compilation  report  and  not  as  a stand-alone  technical  report. 

The  following  component  part  numbers  comprise  the  compilation  report: 

ADPO  13708  thru  ADPO  13761 


UNCLASSIFIED 


Shape-measure  method  for  introducing  the  nearly 

optimal  domain 

A.  Fakharzadeh 

Department  of  Mathematics,  Shahid  Chamran  University  of  Ahvaz . Ahvaz,  Iran. 

a.f  akharzadeh@hotmail . com 

J.  E.  Rubio 

Department  of  Applied  Mathematical  Studies,  University  of  Leeds , Leeds,  LS2  9JT,  UK. 

Abstract 

We  deal  with  introducing  a new  algorithm  for  solving  the  optimal  shape  problems  in 
which  they  are  defined  with  respect  to  a pair  of  geometrical  elements.  The  problem  is 
to  find  the  optimal  domain  approximately  for  a given  functional  that  is  involved  with 
the  solution  of  a linear  or  nonlinear  elliptic  equation  with  a boundary  condition  over 
a domain.  The  Shape-Measure  method,  in  Cartesian  coordinates,  will  be  used  to  find 
the  nearly  optimal  solution  in  two  steps.  By  transferring  the  problem  into  a measure- 
theoretical  form,  first  we  will  find  the  solution  of  the  elliptic  problem  for  a given  domain 
by  using  the  embedding  method.  Then  the  Shape-Measure  method  will  be  applied  to  find 
the  best  domain  approximately.  An  example  will  be  given. 


1 Introduction  and  Problem 

Consider  the  optimal  shape  (optimal  shape  design)  problems  in  which  they  are  defined 
with  respect  to  a pair  of  geometrical  elements;  this  pair  consists  of  a measurable  set 
(in  IR2),  which  can  be  regarded  as  a domain,  and  a simple  closed  curve  containing 
a given  point,  which  is  the  boundary  of  the  set.  By  considering  the  property  for  the 
desired  curves  to  be  simple,  the  problem  depends  on  the  geometry  which  is  used.  In 
polar  coordinates,  we  solved  the  similar  problem  in  [1],  But  in  Cartesian  coordinates, 
it  is  difficult  to  introduce  a linear  condition  which  determines  the  property  of  a closed 
curve  being  simple.  Thus  here  we  consider  some  limitation  on  shape  in  order  to  make 
sure  that  it  is  simple.  The  problem  will  be  solved  in  two  stages.  First,  by  use  of  measures, 
the  value  of  the  objective  function  will  be  calculated  for  any  given  domain.  Then  the 
optimal  domain  will  be  obtained  by  use  of  optimization  techniques. 

Let  D C IR2  be  a bounded  domain  with  a piecewise-smooth,  closed  and  simple 
boundary  dD.  We  assume  that  some  part  of  d D is  fixed  and  the  rest,  IT,  with  the  given 
initial  and  final  points  A and  B respectively,  is  not  fixed.  Here  we  suppose  that  the  fixed 
part  of  3D  is  made  by  three  segments,  parts  of  lines  y = 0.  x = 0 and  y = 1 between 
points  .4(1,0),  (0, 0),  (0, 1),  B(l,  1)  (see  Figure  1). 

Thus  T is  chosen  as  an  appropriate  variable  curve  joining  A and  B so  that  D is 
well-defined.  Let  u( X)  (A  = ( x,y ) G IR2)  be  a bounded  solution  of  the  following  elliptic 
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Fig.  1.  An  admissible  domain  D under  the  assumptions  of  the  numerical  work, 
equation: 

Au(X)  + f(X,u)  = v(X),uleD=0,  (1.1) 

where  XeD-t  v(X)  € JR  is  a bounded  real  function  (v  also  can  be  considered  as  a 
fixed  control  function);  the  function  / is  assumed  to  be  a bounded  and  continuous  real- 
valued function  in  L2(D  x JR).  Moreover  the  above  domain  D is  called  an  admissible  if 
the  equation  (1.1)  has  a bounded  solution  on  D;  we  denote  by  V as  the  set  of  all  such 
admissible  domains.  We  are  going  to  solve  the  problem  of  minimizing  the  functional 
1(D)  = fD  f0(X,u)  dX,  on  the  set  V where  f0  is  a given  continuous,  nonnegative,  real- 
valued function  on  D x JR.  To  calculate  the  value  of  1(D)  for  a given  domain  D,  it  is 
necessary  first  to  identify  the  solution  of  (1.1). 

2 Weak  solution  and  metamorphosis 

In  general,  it  is  difficult  to  identify  a classical  solution  for  the  problem  like  (1.1)  and 
usually  one  tries  to  find  a weak  (generalized)  solution  of  them.  Hence  the  variational 
form  of  (1.1)  is  introduced  in  the  following;  we  remind  the  reader  that  Hq(D)  = 
{ip  e Hl(D)  : ip\BD  = 0} , where  H1(D)  is  the  Sobolev  space  of  order  1. 

Proposition  2.1  Let  u be  the  classical  solution  of  (1.1),  then  we  have  the  following 
equality: 

f (uAip  + iff)  dX=  [ ipv  dX  , 'iip  G Hl(D).  (2.1) 

Jd  Jd 

Proof:  Multiplying  (1.1)  by  the  function  ip  e Dq(D)  and  then  integrating  over  D,  with 
use  of  the  Green’s  formula  (see  [3])  gives  JD(uAip  + ipf-  ipv)  dX  = f9D(ip§%  - u|^)  dS, 
where  n is  the  unit  vector  normal  to  the  boundary  8D  and  directed  outward  with  respect 
to  D.  Because  ip\gD  = 0 and  U|SD  = 0,  then  (2.1)  is  satisfied.  □ 
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Definition  2.2  A function  u £ Hl  (D)  is  called  a bounded  weak  solution  of  the  problem 
(1.1)  when  it  is  bounded  and  satisfies  the  equality  (2.1)  for  all'll)  £ Dq(D)  (the  conditions 
for  existence  of  the  weak  solution  of  the  problem  (1.1)  and  also  the  boundedness  property 
of  it,  have  been  considered  in  many  references,  like  [3]  and  [2]). 

Now  we  apply  our  new  way  which  is  called  the  Shape-Measure  method.  Let  Q = UxD, 
where  U C 1R  is  the  smallest  bounded  set  in  which  the  bounded  weak  solution  «(•) 
takes  values.  Then  by  applying  the  Riesz  Representation  Theorem  ([6]),  a bounded 
weak  solution  can  be  represented  by  a positive  Radon  measure;  the  proof  of  the  following 
Proposition  is  similar  to  the  Proposition  3.1  in  [1], 

Proposition  2.3  Let  u(X)  be  a bounded  generalized  solution  of  (1.1);  there  exist  a 
unique  positive  Radon  measure,  say  pu,  in  A4+(fl)  such  that: 

pu(F)  = f Fdpu  = [ F(X,u)dX-  VFeC(fi).  (2.2) 

Jn  J d 

Thus  the  equality  (2.1)  can  be  changed  to  pu(Fy,)  =7 4,,  V'tp  £ H) (D),  where  Fy,  — 
uAip  + ftp  and  7y,  = fD  ifv  dX.  Also,  1(D)  = pu(f0).  Because  the  measure  pu  projects 
on  the  (x,y)- space  as  the  respective  Lebesgue  measure,  we  should  have  pu(£)  = fle, 
where  £ : — ► 1R  depends  only  on  variable  X (i.e.  £ £ Cj(fl)),  and  Of  is  the  Lebesgue 

integral  of  £ over  D.  Therefore  the  original  problem  can  be  described  as  follows: 

To  find  a measure  pu  € A4+(fi)  so  that  it  satisfies  the  following  constraints: 

M-FV)  = W>e#o(D); 

MO  = V£eCMfi).  (2.3) 

As  Rubio  did  in  [5],  to  be  sure  that  we  do  not  miss  any  solution,  we  extend  the  underlying 
space;  instead  of  finding  a measure  pv  £ introduced  by  Proposition  2.3  and 

equalities  (2.3),  we  seek  a measure  p £ M+(Ll)  which  satisfies  just  the  conditions: 

f*(Fil>)  = £ Hq(D); 

= as,  V£GC,(S1).  (2.4) 

3 Approximation 

The  system  (2.4)  is  linear  because  all  the  functions  in  the  right-hand-side  of  equations 
are  linear  functions  in  their  argument  p.  But  the  number  of  equations  and  the  underlying 
space  are  not  finite.  We  shall  develop  this  system  by  requiring  that  only  a finite  number 
of  the  constraints  are  satisfied.  This  will  be  achieved  by  choosing  countable  sets  of 
functions  whose  linear  combinations  are  dense  in  the  appropriate  spaces.  But  first  we 
should  approximate  the  unknown  part  of  the  boundary  just  by  the  finite  number  of  its 
points.  This  idea  comes  from  the  approximation  of  a curve  by  broken  lines.  For  the  given 
D and  hence  for  the  given  T,  let  Am  = (xrn , ym),  m = 0,1.2, ... , M.  be  a finite  number 
of  points  on  T (where  A0  = A).  We  link  together  each  pair  of  consecutive  points  Am 
and  Am+ 1 for  m = 0, 1, . . . , M — 1 and  close  this  curve  by  joining  the  points  Am  and 
B together.  Now  the  resulted  shape,  which  is  denoted  by  8Dm,  is  an  approximation  for 
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dD;  we  also  call  the  domain  which  introduced  by  its  boundary  8Dm  as  Dm  (see  Figure 
1). 

It  is  possible  that  by  increasing  M,  the  curve  ODm  will  become  closer  and  closer  (in 
the  Euclidean  metric)  to  the  curve  dD,  and  hence  one  may  conclude  that  the  minimizer  of 
I over  T>m  , if  it  exists,  tends  to  the  minimizer  of  I over  T>,  if  it  exists.  But  some  difficulties 
could  arise  (too  oscillatory  a curve  may  cause  problems).  Thus,  we  will  fix  the  number  of 
points.  For  a given  M,  let  the  value  of  the  components  J/i , J/2 >•••■>  Vm,  be  fixed.  Because 
xm  is  a free  term,  the  point  Am  could  be  anywhere  on  the  line  y = Ym,x  > 0 for  every 
m (see  Figure  1).  Therefore  points  Am  and  Am+ 1 can  be  chosen  so  that  they  belong  to 
r and  hence  the  part  of  T between  the  lines  y = Yrn  and  y = Ym+i  can  be  approximated 
by  the  segment  AmAm+ 1.  Hence,  we  do  not  lose  generality.  Thus,  we  fix  the  components 
yi,  j/2)  • • • , Vm  with  the  values  Yj,  Y2, . ■ . , Ym,  respectively. 

Now  we  introduce  the  set  {?/>*€  Hq{D)  : i — 1,2, ... } so  that  the  linear  combin- 
ations of  the  functions  {ipi}  are  uniformly  dense  (that  is,  dense  in  the  topology  of 
the  uniform  convergence)  in  Hq(D).  We  know  that  the  vector  space  of  polynomials 
with  the  variable  x and  y,  P(x,y),  is  dense  in  C°°(D);  therefore  the  set  Po(x,  y)  = 
{p(x,y)  £ P(x,y)Jj>(x,y)  = 0,V(x,t/)  £ dD}  , is  dense  (uniformly)  in  {h  eC°°(D)  : 
h|SD  = 0}  = Cq°(jD)}.  Since  the  set 

Q(x,y)  = {l,x,y,x2,xy,y2,x3,x2y,xy2,y3,...} 

is  a countable  base  for  the  vector  space  P(x,  y),  each  elements  of  P(x,  y)  and  also  Po{x,  y), 
is  a linear  combination  of  the  elements  in  Q(x,y).  By  Theorem  3 of  Mikhailov  [3]  page 
131,  the  space  C°°{D)  is  dense  in  Hl(D);  thus  the  space  C§°(D)  will  be  dense  in  Hq(D). 
Consequently,  the  space  P0(x,y)  is  uniformly  dense  in  Hq(D).  We  define 

M 

tpi(x,y)  = xy(y-l)~[[(x-xi  + y-Yi)qi(x,y),  (3.1) 

(=1 

where  £ Q(x,  y).  Therefore  ^|r  = 0 and  the  set  {ijy(x,  y)  : i = 1, 2, . . . } , is  total  (dense 
in  the  topology  of  the  uniform  convergence)  in  Hq(D). 

For  the  second  set  of  functions,  let  L be  a given  positive  integer  and  divide  D into  L 
(not  necessary  equal)  parts  D\,  D2,  ■ ■ ■ , Dl,  so  that  by  increasing  L the  area  of  Ds,s  = 
1,2 ,L,  will  be  decreased.  Then,  for  each  s we  define: 


£s{x,y,u 


1 

0 


if  ( x,y ) G Ds, 
otherwise. 


These  functions  are  not  continuous,  but  each  of  them  is  the  limit  of  an  increasing  sequence 
of  positive  continuous  functions,  {£•!*};  then  if  p is  any  positive  Radon  measure  on  0, 
/x(Cs)  = limfc-Hx,  p(£Sk)-  The  linear  combination  of  functions  : j = 1, 2, . . . ,L}  for  all 
positive  integer  L,  can  approximate  a function  in  C\  (0)  arbitrary  well  (see  [5]  Chapter 
5). 

By  selecting  just  the  finite  number  of  functions  in  the  mentioned  spaces  the  problem 
(2.4)  can  be  replaced  by  another  one  in  which  we  are  looking  for  the  measure  pMltM2  G 
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A4+(fi),  so  that  it  satisfies  the  following  constraints: 

FM1MAFi)  = 7*.  i = 1,2,...,  Mi; 

MMi,a/2(0)  = aF  j — 1)2,...,  M2,  (3.2) 

where  Mi  and  M2  are  two  positive  integers  and  F = , 7*  = 7^,4  , aj  = a$r  If 

we  denote  by  Q(Mi,M2)  the  set  of  positive  Radon  measures  in  A4+(fi)  which  satisfy 
equalities  (3.2),  and  also  denote  by  Q the  set  of  positive  Radon  measures  in  A4+(fl)  which 
satisfy  equalities  (2.4),  one  can  easily  prove  the  following  Proposition  by  considering  the 
proof  of  Proposition  I II.  I in  [5]. 

Proposition  3.1  : If  Mj,  M2  — > oo  then  Q(Mi,M2 ) — > Q;  hence  for  the  large 
enough  numbers  Mi  and  M2  the  set  Q can  be  identified  by  Q(Mi,M2). 

But  even  if  the  number  of  equations  in  (3.2)  is  finite,  the  underlying  space  Q{M\,  M2) 
is  still  infinite-dimensional.  By  Theorem  A. 5 in  the  Appendix  of  [5],  tLMl,M2  in  (3.2) 
can  be  characterized  as  g,MuM2  = ocnS(Zn),  with  triples  Z„  £ Q and  the 

coefficients  an  > 0 for  n — 1,2, ...,M\  + M2,  where  S(z)  £ A4+(fi)  is  supposed  to 
be  a unitary  atomic  measure  with  support  the  singleton  set  {z}.  Thus  the  measure 
problem  is  equivalent  to  a nonlinear  one  in  which  the  unknowns  are  the  coefficients  an 
and  supports  {Znj.  Proposition  III. 3 of  [5]  Chapter  3,  states  that  the  measure 
has  the  following  form 

N 

Fmum2  ~ ^ ) On d(ZTI ) , (*I-^) 

71=  1 

where  Zn,  n = 1.2, ... , N,  belongs  to  a dense  subset  of  fi.  Now  let  us  put  a discretization 
on  fi,  with  the  nodes  Z„  = (x„,  yn  ,un),  in  a dense  subset  of  ft;  then  we  can  set  up  the 
following  linear  system  in  which  the  unknowns  are  the  coefficients  an: 

an  > 0,  n = 1, 2, . . . , N\ 

N 

^ OinFifZjf)  ~7 ii  i = 1, 2, . . . , Mi , 

71  = 1 

N 

^ an£j(Zn)  — flj)  1 = 1,2,...,  M2.  (fi-1!) 

n—  1 

The  solution  of  (3.4)  is  not  necessary  unique  (even  if  the  problem  (1.1)  satisfies 
the  necessary  conditions  for  having  a unique  bounded  weak  solution),  because  of  the 
approximation  scheme. 

4 The  optimal  solution 

The  main  aim  of  the  present  section  is  to  find  an  optimal  domain  D*  £ Dm  so  that 
the  value  of  I(£>*)  will  be  the  minimum  on  the  set  VM.  By  applying  the  result  of  the 
previous  section,  a solution  of  (1.1)  can  be  found.  Indeed,  it  is  approximated  by  a solution 
of  the  linear  system  (3.4)  according  to  the  variables,  xm,  m = 1, 2, . . . , M.  As  mentioned, 
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this  solution  is  not  necessary  unique.  Let  us  specify  one  by  solving  the  following  linear 
programming  problem 

N 

n= 1 

an  > 0,  n—  1,2,..., AT; 

N 

^ ' ^n-S'i^n)  = T i)  % ~ 1,  2,  . . . , M\] 

n= 1 
iV 

^ ' &n£j(Zn)  — a.j,  j = 1, 2, . . . , ili^.  (4-1) 

n=l 

Thus,  for  each  D,  the  value  1(D)  = JDf0(X,u ) dX  = p{f0)  — /*m,,m2(/o),  is  defined 
uniquely  in  terms  of  the  variables  xm,  m = 1, 2, . . . , M.  So,  we  set  up  a function,  J,  on 
VM  defined  by  D £ VM  — > 1(D)  = fiMuM2(fo)  where  Mm,,m2(/o)  = E!Li  atnf0(Zn). 
Clearly  J can  be  regarded  as  a vector  function: 

J : (xi,x2,...,xM)  € JRm  — *•  (J,MuM2(fo)  € 1R.  (4.2) 

Since  J is  a real-valued  function  which  is  bounded  below,  and  is  defined  on  a compact 
set  (since  constraints  are  to  be  put  in  the  variables),  it  is  possible  to  find  a sequence  of 
points  so  that  the  value  of  the  function  along  the  sequence  tends  to  the  (finite)  infimnm 
of  the  function.  The  coordinate  values  corresponding  to  the  points  in  the  sequence  are  of 
course  finite.  Now,  suppose  that  (x\,  X2i  ■ ■ 'i*m)  is  the  minimizer  of  the  vector  function 
J;  it  can  be  identified  by  using  one  of  the  related  minimization  methods.  The  introduced 
domain  by  the  minimizer  is  denoted  by  D* . We  assume  in  the  following  theoretical  result 
that  the  minimization  algorithm  which  is  used,  is  perfect;  that  is,  it  comes  out  with  the 
global  minimum  of  J in  its  (compact)  domain. 


Theorem  4.1  : Let  M,  M\  and  M2  be  the  given  positive  integers  which  were  defined  in 
section  3,  and  D*  be  the  minimizer  of  (4-2)  as  mentioned  above.  Then  D*  is  the  min- 
imizer domain  of  the  functional  I over  T>m  and  the  value  o/I(D*)  can  be  approximated 
by  J(D*);  moreover  J (D*)  — » I(D*)  as  and  M2  tend  to  infinity. 

Proof:  Suppose  D*  is  not  the  minimizer  of  I;  hence  there  exists  a domain,  call  D', 
in  T>m  so  that  I(D')  < I(D*).  Proposition  2.3  shows  that  there  is  a unique  measure, 
call  p1,  in  A4+(D)  so  that  I(D')  = fi'(fo),  and  also  Proposition  3.1  states  that  for 
sufficiently  large  numbers  Mi  and  M2,  fx’(Jo)  can  be  approximated  by  p'Mi  M2(f0 ) in 
Q(Mi,M2).  Thus,  I(D')  = Mmi ,m2(/°)  = In  the  same  way,  one  can  show  that 

J(D*)  approximates  I(D*);  so  I(D*)  = p*Mi  M2(fo)  = J(D*).  Hence  J(D')  < J(D*). 
which  is  contrary  with  the  fact  that  D*  is  the  minimizer  of  J.  Moreover,  from  Proposition 
3.1  it  follows  that  J(D*)  tends  to  I(D*)  as  Mj, M2  — * 00.  □ 
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5 Numerical  example 

We  consider  the  elliptic  equations  (1.1)  with 

V(XV)  = { 1 if  (x,y)eDnC, 

' ’ \ 0 otherwise, 

where  C is  the  square  [|,  |]  x [|,  |]  (see  Figure  1).  We  also  take  M = 8,  M\  = 10,  M2  = 8, 
N = 740  (the  36  number  of  nodes  are  chosen  so  that  u\0D  = 0)  and  suppose  Y),  Y2,  . • ■ , Fg 
are  0.15, 0.25, . . . , 0.85,  respectively.  By  extra  constraints,  xm  > | , m = 2, 3, . . . , 7,  the 

3 3 

value  of  7 i for  any  D 6 T>m  is  defined  as  7,  = sm  ipi(x,  y ) dxdy  , z = 1,2,...,  10.  We 

4 4 

also  assume  that  the  function  u takes  values  in  U =[-1,1],  and  consider  the  polynomials 
qi(x,  y)  as  l,a :,y,x2,xy,y2,x3,x2y,xy2,y3.  The  function  /0  is  chosen  as  f0  = (u-0.1)2. 
This  function  can  be  considered  as  a distribution  of  heat  in  the  surface  for  the  system 
governed  by  an  elliptic  equations. 

In  minimization,  we  apply  the  Downhill  Simplex  Method  in  Multidimension  by  using 
the  Subroutine  AMOEBA  (see  [4])  and  also  consider  an  upper  bound  for  variables 
(suppose  they  are  not  higher  than  2).  These  conditions  are  applied  by  means  of  the 
penalty  method  (see  [7]).  Hence,  for  the  nonlinear  case  of  the  partial  differential  equations 
(1.1),  we  have  taken  f(x,y,u ) = 0.25w2,  and  used  the  initial  values  as  Xm  = 1.0, m = 
1, 2, . . . , 8,  and  the  stopping  tolerance  for  the  program  (variable  ftol  in  the  Subroutine 
AMOEBA)  as  10~7.  We  remind  the  reader  that  the  functions  F,  and  the  values  of  7 
i = l,2,...',  10,  have  been  calculated  by  the  package  uMaple  V.T . The  results  are:  the 
optimal  value  of  I = 0.45467920356379,  the  number  of  iterations  = 502,  the  value 
of  the  variables  in  the  final  step  are  = 1.05019,  X2  = 1.08521,  X3  = 0.750001, 
X4  = 0.768701,  Xh  = 1.12986,  X6  = 1.13775  ,X7  = 0.97783,  X8  = 1.61566,  which 
represent  the  optimal  domain,  shown  in  the  Figure  2. 
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Fig.  2.  The  initial  and  the  optimal  domain  for  nonlinear  case  of  elliptic  equations. 


